The [1,2] Pade Amplitudes for tttt Scatterings 
in Chiral Perturbation Theory 



GuANG-You QiN, W. Z. Deng, Z. G. Xiao and H. Q. Zheng 
Department of Physics, Peking University, Beijing 100871, P. R. China 

Abstract 

A detailed analysis to the [1,2] Pade approximation to the vrvr scat- 
tering 2-loop amplitudes in chiral perturbation theory is made. 
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The chiral perturbation theory (ChPT) is a powerful tool in 

studying strong interaction physics at low energies and has been extensively 
studied at 1-loop level 0]. The 2-loop results are also available in recent 
years However, since the chiral expansion is an expansion in terms 

of the external momentum, the perturbation series to any finite order di- 
verges very rapidly at high energies. Therefore the violation of unitarity 
gets even worse for 2-loop amplitudes than the 1-loop amplitudes at high 
energies. Also the number of parameters in the effective Lagrangian which 
are not fixed by symmetry alone increases rapidly. Therefore increasing the 
order of the perturbation expansion does not work at all for the purpose of 
exploring physics in the non-perturbative region, or at higher energies and 
non-perturbative studies become necessary. A widely used method to rem- 
edy the violation of unitarity is the so called Pade approximation. [| A nice 
feature of the Pade approximation is that it restores unitarity with full re- 
spect, at low energies, to the available information from perturbation theory. 
Therefore, even though it is well known that it violates crossing symmetry, 
Pade approximation is considered to be a valuable tool in exploring physics 
in the non-perturbative region, such as the properties of physical resonances. 



^Variations of the Pade approximation method can be found in the Hterature named as 
the inverse amphtude method (lAM) and the chiral unitarization approach ||ll| with 
somewhat different formahsm and motivation. 
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However, a previous study |]T2[ indicates that the [1,1] Pade approximation 



encounters a serious problem by predicting spurious physical sheet resonances 
(SPSRs). Usually these SPSRs locate at distant places very far from the re- 
gion where the perturbation results are valid. The predictions of the Pade 
approximants constructed from the perturbative amplitudes should, of course 
not, be considered as meaningful in the region far away from the region where 
perturbation theory remains to be valid. One may further argue that since 
those SPSRs are far from the region we are concerning the use of the Pade 
approximation is still acceptable at least in phenomenological discussions. 
However, the problem with Pade approximation is not only because it pre- 
dicts SPSRs in the distant region too far away to be worthwhile to pay any 
attention, but also because those SPSRs usually have large couplings to tttt 
which lead to strong influence to the region we are really interested in, and 
hence their existence casts doubt on the remaining predictions of the Pade 
approximants which might otherwise be assumed as meaningful, at least at 
quantitative level. The aim of the present study is to further investigate the 
Pade approximation following the method of Ref. [|r^. We will extend the 
work of Ref. [1^ by also analyzing the [1,2] Pade approximants, since the 
2-loop perturbation results are already available. One of the main motiva- 
tion of the present work is to investigate the possibility that the [1,2] Pade 
approximants can rescue, to some extent, the bad situation the [1,1] Pade 
approximants encounter. The conclusion is rather negative, as we will see 
in the following text. However, we believe it is still worthwhile to present 
our results. Since the Pade approximation is a very popular approximation 
method widely used in phenomenological discussions, we hope the presen- 
tation of the present work could benefit physicists who are working in the 
related fields of non-perturbative physics. 

For the tttt Tin scattering, it is well known that the isospin amplitudes 
in the s channel can be decomposed as, 

T^=°(s, t, u) = 3A{s, t, u) + A{t, u, s) + A{u, s, t) , 
T^=^{s,t,u) = A{t,u,s) - A{u,s,t) , 

T^=^{s,t,u) = A{t,u,s) + A{u,s,t) , (1) 
where s,t,u are the usual Mandelstam variables, 

s = {pi+p2Y/M'^ , t = {p,-p,f/Ml , u = {p,-p,f/Ml. (2) 
In SU (2) X SU (2) chiral perturbation theory to two loops |Q, the momentum 
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expansion of the amplitudes amounts to a Taylor series in 

..^|f. (3) 

The amplitudes can be expressed in terms of six parameters bi, 62, &3, &5, 



u? 



y4(s, t, u) = X2{s — 1) 

+xl (bi + b2S + 635^ + biit 

+xl + G^^\s, t) + M) 

+xl (b^s^ + besit - 

+xl lF^'\s) + G^^\s,t) + G^^\s,u)] 

+0{xt) . ' (4) 

The expressions of the functions and and the constants bi 

are rather lengthy and we refer to the original work for the details. 
The partial wave expansion of the isospin amplitudes is written as 

T\s, t, u) = 327T ;^(2 J + l)Pj(cos e)Tj{s) . (5) 
J 

The inverse expression is 

^i(^) = 7^ f' d{cose)Pj{cose)T\s,t,u) , 
cos6 = 1 + 



s-4 ' 

u = 4 - s -t . (6) 
The partial wave amplitude in ChPT expanded to O(p^) is, 

Tj{s)=TUs)+TUs)+TUs). (7) 



In Ref. |]T2[ we have discussed the [1,1] Fade approximants of the partial 



wave amplitudes in 1-loop ChPT. To proceed we now construct the [1,2] 
Pade approximants to the partial wave amplitudes in 2-loop ChPT, 

T^^''\s) = ^ , . (8) 



TUs) TUs) ^ \TUs) 



Perturbation theory satisfies the elastic unitarity relation, 

ImTjis) = p{s)\Tj{s)\\ (9) 

at each order of the perturbation expansion in powers of the quark masses 
and external momentum, i.e., 

ImTj.is) = 0, 

ImTj^M = Pis) {Tj^,is)y , 

ImTj^.is) = 2p{s)Tj^,{s)ReTj^,{s) , (10) 



With these relations it is easy to prove that the [1,2] Pade approximant in 
Eq.(Rp satisfies elastic unitarity: 



Pis)\Tr\s)\' . (11) 

In the following we frequently omit the indices /, J of the T matrix for 
simplicity if it causes no confusion. For any given amplitude satisfying single 
channel unitarity, following the method of Refs. [|1^, , we define two real 
analytic functions F and F as 



F{s) 
F{s) 



1 
2 



Sis) 



2ip{s) 



S{s) 
^S{s) - 



Sis) 



'12) 



It is obvious that F and pF are the analytic continuation of cos(25) and 
sin(25), as the scattering S matrix is equal to exp{2i6} in the physical region. 
According to ||T^, |I4|, we have the following dispersion relations for F and 
F: 



sm{2S) 



pF =p{a + J2 



+ 



2ip{si){s 

1 1 Ijn^^, 

n JL s' — s 



E 



1 



-^2^p{zy)s'{zy){s-zy) 

(13) 



and, 



cos(25) 



a 



+ E 



2(s 



+ E 



- 2S'izy){s - zy) 



1 r Im^,^, 
TV Jl s' — s 



(14) 



4 



where a and a are subtraction constants, Si denotes the possible bound state 
pole positions and Pi denotes the corresponding residues of S; denotes 
either the possible resonance pole positions on the second sheet, which are 
grouped into complex conjugated pairs, or the virtual state pole positions 
when Zj^ is real. The integrals in Eqs. (13) and (14) denote the cut con- 



tributions and one subtraction to each integral is understood, according to 
general physical consideration.^ L = {—oo,0] is the left hand cut {l.h.c). 
The discontinuities on the left in Eqs. (0) and (|Tip satisfy the following 



equations 



ImiF(s) = -2p(s)ImiImRT(s) , 

ImiF(s) = 2ImLRe^jT(s) . (15) 

In order to evaluate the values of ImiRe/jT(s) and Ivul^vhrT^s), we need the 
analytical expressions of ReijT(s) and lm.RT{s) which can be derived from 
the expression of the T^^''^\s) in Eq. (13), 



T2{sfa 
Im^T'-Js = 



Re«T[i'2](.) 



(16) 



where 
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a = r2(s)2-T2(s)ReHT4(s)-r2(s)ReRT6(s) + (Refir4(s))2-(Im^T4(s)) 
h = T2{s)\mRT^{s)+T2{s)\mRT(,{s)-2ReRT^{s)lmRT^{s) . (17) 

Analytical expressions for Im^F and Im^F in terms of perturbation am- 
plitudes can also be written down, or can be calculated numerically from 
Eq. (p!6|). In Eqs. (|1^) and (0) we did not include resonance poles on the 
first sheet, since they are not allowed physically. However, as we stated be- 
fore, the Pade amplitude may contain SPSRs. When using Eqs. ( [TB| ) and 
([T^ ) to analyze the Fade amplitude the dispersion representations have to 
be modified to include those terms representing SPSRs. This can easily be 
done by using Eq. (|1^). 

^In general the dispersion integrals in Eqs. ( [l3| ) and (|lj) need one subtraction, but 
the Pade amplitude is special in that the integrals are finite and need no subtraction. 
Therefore the subtraction constant, a and a are and the integrals are unsubtracted 
when analyzing the Pade amplitude. 
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Making use of the property of the scattering amphtude at threshold one 
can recast Eqs. (|1^) and (|1^ in the following form: 



sm{26)=p{s)F =pis){2a^j + Y: ^'^^ 



2ip{si){s - Si) (4 - Si) 



y 2ip{zy)S'{zl%s - )(4 - ) ^ n Jl (s' - s)(s' - 4) ' 

_ p _ 1 , A (4 - s) ^ 4^s 

cos^^dj - ^ - ' + 4- 2(s - s,)(4 - s,) + ^ 2S'{zj^){s - )(4 - z^^ 



in which Oj represents the scattering length parameter in the corresponding 



channel. The difference between Eqs. (|13|), (|1J) and Eq. ([18| ) really makes 
the difference: in the latter formula a constant contribution is subtracted 
from each pole term. The new definition of the pole contribution (that is the 
original pole contribution minus the its contribution at s = 4) only probes 
the s dependence of the pole term. For example, in the limit oo 
while S'{Zj^)zj^ is held fixed, the pole contributes a constant term to the 
dispersion relation according to Eqs. (0) and This constant term is 

reabsorbed into the scattering length parameter in Eq. (|18D and the pole no 
longer contributes to the dispersion relation, according to the new definition 
of pole contribution. Similar discussion can be made for the case of SPSR. In 
the following we will always use the new definition of the pole contribution. 
Except the pole contributions, the rest of the r.h.s. of Eq. (|18D will be called 
the background contribution in the following text. 

The Eq. (|18D allows us to explicitly examine different contributions from 
various kinds of dynamical singularities to the phase shifts. We have com- 
puted various contributions to cos(2(5) and sin(25) in IJ=00,11 and 20 chan- 
nels both in [1,2] Fade and [1,1] Fade approximations as presented below. 

The SU{2) X SU{2) effective Lagrangian at Oijf') contains two sets of pa- 
rameters: li - U of 0{p'^) and ri - rg of 0{p^). Here we take these parameters 
the same as in Ref. [Q]: the scale-independent couplings k are, 

h = -1.7 , I2 = 6.1 , [3 = 2.9 , U = 4.3 ; (19) 

and the constants (the resonance contributions to the low-energy constants 
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648MeV(M) 


118MeV(r) 
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SPSR 


69.25 


19.22 


-0.14+17.89i 


20 


vs 


0.0513651 





0.0477297 


SPSR 


-18.58 


14.15 


8.44+6. 13i 


SPSR 


135.14 


32.16 


0.38+32.34i 



Table 1: Resonances (R), spurious physical sheet resonances (SPSR) and vir- 
tual states (VS) as predicted by the [1,2] Pade approximation on the com- 
plex s plane using the coupling constants given by Eq. (|TP|) and Eq. (PU]). 
The pole position Sp = (M + ?r/2)^. All numbers are in unit of un- 
less otherwise stated. The values of the scattering length parameters are, 
a° = 0.224, a2 = -0.0412. 



of 0(p^)) are. 



rf = -0.6 X 10~^ 
rf = -1.0 X 10"^ 



rf = 1.3xl0"\ rf = -1.7x10" 
= 1.1 X 10"^ , rf = 0.3 x 10^^ 



(20) 



and we take the renormalization scale /i = IGeV when evaluating the con- 
stants hi appeared in Eq. (^). Using the above values of parameters we can 
determine poles and cuts of the Pade amplitudes. As shown in Table |I], Pade 
approximation not only predicts the existence of the a and p resonances, 
but also generates many other poles on the complex s plane, and more poles 
exist in [1,2] than in [1,1] Pade approximant. In addition to the poles found 
in table |l|, it is found that there also exist 2 pairs of BS/VS poles located 
close to the Adler Zero position of T^^ iii both IJ=00 and IJ=20 channels 
(in a wide range of the li and rj parameters). Similar to what happens in the 
[1,1] Pade case [0] they can be tuned away within reasonable range of the 
li and Ti parameters and are only artifacts of the Pade approximants. More 
importantly they only have very tiny effects and can be safely neglected. 
The existence of the virtual state in the IJ=20 channel has been clarified in 
Ref. [0 but its effect is also very small. 

Besides those well established resonances which can be found in table 1, 
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Fade approximants predict resonances or SPSRs at distant places on the 
complex s plane. Chiral perturbation expansion only works in a region close 
to the threshold or \s\ << IGeV^, therefore the Fade approximants should 
also be expected to be reasonable only in a limited region: \z\ << IGeV^ 
on the complex s plane. Hence any prediction from Fade approximants at 
distant places should not be trustworthy, no matter the predicted poles are 
on the first sheet or on the second sheet. Of course we should not take 
these predictions seriously. The real problem for the Fade amplitude is that 
in many cases the distant poles do not truly decouple from the low energy 
physics, as indicated by their large couplings. The Eqs. ([T3|) and ( |l4|) afford 



us a useful tool to evaluate the influence of these spurious or unreliable 
contributions quantitatively. 

To have a clear insight to the problem we are facing we perform the 
following calculation: We use the MINUIT program to make a global fit to 
the experimental phase shift in both the IJ=00,20 and 11 channels using 
the Fade amplitudes. The data ([[15|-|T^) are taken from the threshold to 
730MeV in IJ=00 channel and to IGeV in IJ=11 and 20 channels. The fit 
results of the coupling constants k and are listed in the following, 



h = 


—4 47+0-82 

^•^'-0.83 ' 


l2 = 


4 Q7+0.11 
^•'J ' -0.11 




h = 


r, 91 +14.00 
'-'•^-'--13.71 5 


k = 


7 qc+0.51 
' •'-'^-0.52 




ri = 


-84.851^1:1 


X 10" 




3Mtllf, X 10-^ 


^^3 = 


-34.3611^^ X 


10"^ 


ri = 


-o.utlit X 10-^ 


r5 = 


7.661}:?^ X 10 


-4 

5 


re = 


-2.6lt°:i X 10-^ 



(21) 

From above results we find that the fit is very insensitive to I3, ri and 
r2. The other parameters are barely comparable in order of magnitude 
to those given in Eq. ( PP] ) p(]| . The results for the fit are shown in table 



and Figs. [l|-0. In order to compare the [1,2] Fade results with [1,1] Fade 
results we also made the global fit in the latter case. The results are given 
in table 0. In Fig. ^ we plot the [1,1] Fade fit results in the most interest- 
ing IJ=00 channel. In Figs. 1-4, the solid lines represent the phase shift 
value directly obtained from the Fade amplitudes in the physical region, i.e., 
= Atan(ImRT(s)/Rei?T(s)). The dash-dot-dot-dash lines in Figs. 1-3 
and the dash-dot-dash line in Fig. 4 represent the summation of all contribu- 
tions on the r.h.s. of Eqs. (|1^) and (|1^) (also including SFSRs' contributions). 
The two kinds of lines must coincide with each other as a consistency check 
of our numerical calculation. The coupling constants obtained from the fit 
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Table 2: Resonances (R), spurious physical sheet resonances (SPSR) and vir- 
tual states (VS) as predicted by the [1,2] Pade approximation on the complex 
s plane using the values from Eq. (|21|). The pole position Sp = (M + ?r/2)^. 
All numbers in the table are in unit of unless stated otherwise The fit 
values of the scattering length parameters are, = 0.217, = —0.0531. 

using the [1,1] Pade amplitudes are the following, 

h = — 0.33^o;Qgg , I2 = 5.83^o;Qg7 , I3 = 24.09lj;|!,;52 , h = 3.51^o;27 . (22) 

Comparing the results of [1,1] Pade approximants and the [1,2] Pade ap- 
proximants we conclude that the [1,2] Pade approximation does not in general 
improve the bad situation the [1,1] Pade approximants encounter, though the 
former can give a better fit to the phase shift data^. In the IJ=11 channel 
the [1,1] Pade amplitude only predict the p resonance, in the [1,2] amplitude 
there are additional poles at distant positions but their contributions are 
rather small, so both the two amplitudes can be considered as phenomeno- 
logically successful. In the IJ=20 channel both amplitudes run into disaster 
since both of them predict huge contributions from SPSRs. Comparing Fig. g 
with Fig. 2 of Ref. ||12[ one finds that the [1,2] Pade amplitude is even worse in 
the sense that it predicts large SPSR contribution also to sin(25). In fact the 
perturbation result in this channel is much better |^. In the most interesting 
IJ=00 channel, the situation is more complicated. From Fig. ^ we find that 

■^This is hardly surprising since it contains 6 more parameters. 
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0.0318398 
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Table 3: Resonances (R), spurious physical sheet resonances (SPSR) and vir- 
tual states (VS) as predicted by the [1,1] Pade approximation on the com- 
plex s plane using the values from Eq. (|T9D and Eg. (|20|) . The pole position 
Spoie = {M + ir/2)^. All numbers in the table are in unit of unless 
stated otherwise. The fit values of the scattering length parameters are, 
= 0.186, a2 = -0.0467. 

the SPSR's contribution to sin(2(5) is sizable yet in Fig. |]the SPSR's contri- 
bution is very small. However, in the latter case there appears another large 
contribution from a resonance (R) located at sr = (0.318 ± 0.468i)^GeV^ 
in the IJ=00 channel, which is neither very far away from nor very close 
to the low energy physics region. When looking at table one may even 
confuse the two resonances in the IJ=00 channel: the a and R. The two 
resonances are distinguishable in the following way: in table ^ one resonance 
pole is much closer to the physical region comparing with another one, the 
former is of course denoted as a. Q Then we can tune the parameters li and 
Ti continuously, from the values given in Eqs. (|19D and (pO|) to the central 
values of the parameters given in Eq. (^), and keep track of the pole po- 
sitions. In this way we can distinguish the two resonance poles in table |^. 
Obviously the pole position of R is very sensitive to the parameters of the 
chiral Lagrangian and the pole position of the a meson is rather stable. The 
former comes actually from very distant places, therefore its existence and 
contribution is still doubtful even though it is located not very far from the 
physical region as predicted by the parameter set determined from the global 
fit. In this sense one hesitates to conclude that the [1,2] Pade approximant 
in the IJ=00 channel is better than the [1,1] approximation even though the 
SPSRs' contribution is reduced. 

^The question whether this resonance is the a meson responsible for spontaneous chiral 
symmetry breaking is another question, for recent discussions, see for example Refs. 
and 



10 



In above discussions we point out and study in details the problem of 
the Pade approximation method induced by the existence of SPSRs, which 
has not drawn much attention in the literature. It is worth mentioning 
that the author of Ref. |^ noticed such a problem and proposed a scheme 
to remedy the situation when studying the tt scalar form-factor, F, to two 
loops using IAM.0 When obtaining Eq. (7) from Eq. (6) in that paper, the 
author essentially did the following: expanding the SPSR term up to O(s^) 
and neglecting the higher order terms, and reabsorb the coefficients of the 
second order polynomial into the low energy constants, and the latter are 
determined by either CHPT or by fit. In this way the SPSR (which, according 
to the author, locates on the negative real axis) is eliminated in the newly 
obtained F. This effort is respectable as it is in the right direction, but is still 
problematic. Actually it is not difficult to prove that the procedure proposed 
by the author does not make the effects of the spurious pole totally vanish 
but actually moves the pole from the 1st sheet to the 2nd sheet (see the 
appendix for the proof ).Q A pole on the negative real axis on the 2nd sheet 
is still dubious, even though its numerical influence is not estimated. 

To conclude, in general the prediction from the [1,2] Pade approximation 
is not in any sense more trustworthy than the [1,1] Pade approximation. A 
lesson one may draw from the discussion made in this paper is that physics at 
distant places (no matter spurious or not) as predicted by the Pade approxi- 
mation do not necessarily decouple at low energies. We suggest, to make safe 
use of the Pade approximation one has to make a case by case analysis to the 
amplitudes using the method proposed in this paper. The smallness of the 
contribution from high energies may be considered as a necessary condition 
for the predictions of the Pade amplitude at moderately low energies to be 
numerically trustworthy. 

This work is supported in part by China National Natural Science Foun- 
dation under grant number 10047003 and 10055003. 
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A Appendix 



The proof follows: denote the formfactor defined in Eq. (6) in Ref. by 



Fi and that defined in Eq. (7) by F. The major difference between F and 
Fi is that the latter contains the SPSR. Especially we have ImF=ImFi and 

F = P^is) + - / ds'——^^ = ReF{s) + ilmF{s) . (23) 

This F, even though does not satisfy exact unitarity, maintains the two-sheet 
structure and can be analytically continued to the second sheet by changing 



the sign of the kinematic factor p = y 1 — Am'^/s, notice that both Re-F(s) 
and ilmF{s) can be analytically continued to the complex s plane. To be 
specific F^^ is, 

F'' = pJs) + -P r ds' -ilmFis) = ReFis)-i\mFis) , (24) 

vr J Ami S''^[S' — S) 

i.e., the F^^ is obtainable by just simply changing the sign of lmF{s) and 
P stands for the principle value. By construction F does not contain any 
poles on the entire cut plane and according to Ref. |^ F^^ does. That is 
possible only when both ReF(s) and ilmF{s) contain poles on the cut plane 
but cancel each other in ReF(s) + ilmF{s), but not in ReF(s) — ilm.F{s). 
Therefore the pole locations of F^^ are the same as ilmF{s) and are also the 
same as ilmFi{s). However it is clear from Eq. (6) of Ref. p3[ that ilmFi{s) 
contains the poles of Fi on the 1st sheet (corresponding to zeros of T^'^\ for 
the definition of the latter see Ref. [^), as well as poles on the 2nd sheet 
(corresponding to zeros of F*^^)*). Therefore all poles of Fi are transmitted 
into F^^ including both the a pole and the spurious pole. 
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Figure 1: Various contributions to cos(25q) and sin(2(5Q) in the IJ=00 chan- 
nel using the values from Eq.fpTl). Notice that in the IJ=00 channel we only 
fit the data up to 730MeV. The [1,2] Fade amplitude gives a much better 
description to the data above 730MeV than the [1,1] Fade amplitude, see 
Fig. ^ for comparison. The solid line and the dash-dot-dot-dash line are ex- 
plained in the text. The dashed line is from the background contribution 
only, whereas the dotted line represents the contributions from the a and the 
background, and the dash-dot-dash line is to add the resonance (R) contri- 
bution to the dotted line. The fact that the dash-dot-dash line is very close 
to the solid line indicates that the SFSR's contribution is very small. 




Figure 2: Various contributions to cos(25j) and sin(2(5j) in the IJ=11 chan- 
nel using the values from Eq.(pID. The meaning of different lines is similar 
to those in Fig. 1. 
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Figure 3: Various contributions to cos(2(5q) and sin(25Q) in the IJ=20 chan- 
nel using the values from Eq. (0). The virtual state contribution as repre- 
sented by the difference between the dashed line and the dotted line is very 
small. 



0.5 



O 

O -0.5 J 



-1.0 




0.2 



o 

to 



1-5-1 



1.0- 



0.5- 



0.0- 



-0.5- 



-1.0- 



0.2 



- Pade 
backgroud contribution 
-with a resonance 
■with all poles 



0.3 



0.4 



0.5 

s'''(GeV) 



-Pade'^'^' 

backgroud contribution 
■ with a resonance 
■with all poles 



0.3 



0.4 



0.5 

s'''(GeV) 



0.6 



0.7 




0.6 



0.7 



Figure 4: Various contributions to cos(25q) and sin(25Q) in the IJ=00 chan- 
nel using the values from Eq. (^2]), using the [1,1] Fade amplitude. 



